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Summary 

Recently developed traction data for Rocket Propellant 1 
(RP-l), a hydrocarbon fuel of the kerosene family, were used 
to develop the parameters needed by the bearing code 
SHABERTH in order to include RP-l as a lubricant choice. 
As an aid to future additions, a review of the procedure tor 
inputting data for a new lubricant choice is presented. In 
particular, the fluid traction model used by SHABERTH is 
discussed, and issues concerning it are presented. In the 
process of fitting the RP-l traction data to the model, certain 
assumptions and simplifications were necessary. The error 
resulting from making these simplifications is discussed. 

A slight error was detected in the traction equations that are 
in the SHABERTH program, but a lack of traction data on 
the original lubricants prevented immediate correction of the 
error. However, the error was judged insignificant to the 
traction prediction for RP-l . A temperature dependency that 
is not accounted for by the model was found in a particular 
traction coefficient for RP-l. The maximum traction coeffi- 
cients that were predicted by using temperature-averaged 
coefficients erred by an average of 12 percent, whereas those 
predicted by using the discrete coefficients at the two 
temperature conditions erred by an average of 4 percent. 
Currently, the temperature-averaged coefficients are being 
used. The model of the pressure-viscosity coefficient as a 
function of temperature was modified tor RP-l by using 
pressure-viscosity data valid tor RP-l. 


Background 

The introduction of reusable engines tor space launch 
vehicles has made necessary the development of more rigorous 
design methodologies for long-life bearings in launch vehicle 
turbopump applications. Considerations of weight and 
simplicity require that the working fuel and oxidizer be used 
as the turbopump coolant and lubricant; however, there is 
currently a lack of data on the rheological properties of the 
types of fluids used in these applications at the conditions 
present in a bearing contact (i.c., high pressures and potentially 
high shear rates). 

One fuel that has been considered for use in future launch 
vehicles (e.g., the advanced launch system (ALS)) is Rocket 
Propellant 1 (RP-l), a type of kerosene. Though RP-l has been 


used in the past to cool and lubricate the bearings of expendable 
launch vehicles, such as Titan, only a small amount of data 
on its rheological properties has been collected (rets. 1 to 3). 

In order to design long-life bearings to be run in RP-l, the 
bearing behavior needs to be theoretically analyzed. One of 
the most widely used rolling-element-bearing design and 
analysis tools is the computer code SHABERTH (Shaft Bearing 
Thermal Analysis, ref. 4). SHABERTH simulates the 
thermomechanical performance of a load support system 
consisting of up to five ball, cylindrical, or tapered-roller 
bearings. Transient or steady-state temperatures can be 
calculated by using a lumped-mass thermal model that takes 
into account free convection, forced convection, conduction, 
radiation, and mass transport heat transfer. 

Since the SHABERTH code takes into account the influence 
of the lubricating fluid on bearing behavior, it requires certain 
empirical rheological data for the particular fluid. The 
SHABERTH code currently contains hard-coded rheological 
data on a number of lubricants. These data consist of density , 
thermal coefficient of expansion; thermal conductivity; the 
relationships of viscosity and temperature, and piessure- 
viscosity coefficient and temperature at ambient pressure; and 
the characteristics of fluid traction as a function of shear rate. 

The rheological properties of RP-l. a low-viscosity fuel, 
differ greatly from those of the lubricating oils currently in 
the program as lubricant choices; therefore, the SHABER1 H 
code had to be modified to include RP-l . The density, thermal 
expansion, and thermal conductivity, as well as data for 
viscosity with respect to temperature and for the pressure- 
viscosity coefficient as a function of temperature were avail- 
able. The characteristics of traction force as a function of shear 
rate, however, needed to be established. Once all of the 
rheological data were established, they had to be lit to the 
various rheological models used by the code; this determined 
the inherent parameters that needed to be hard-coded into the 
SHABERTH code. 


Introduction 

The purposes of this report are (1) to review the general 
procedure for inputting data for a new' lubricant choice into 
the SHABERTH bearing code and to describe the theoretical 
traction model used by SHABERTH, (2) to describe the luting 
of the experimental RP-l fluid traction data to the model, and 
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(3) to present results and quantify the error due to making 
certain assumptions. AH RP-1 parameters needed by 
SHABERTH to run the program are specified. 

As previously mentioned, in order to predict the behavior 
of a bearing being lubricated and cooled by a specific fluid, 
the SHABERTH bearing code requires empirically derived 
rheological data on the fluid (viscosity-temperature 
relationship, density, thermal expansion, thermal conductivity, 
pressure-viscosity coefficient, and traction force as a function 
of shear rate). The experimental data must be fit to various 
rheological models that the code uses to determine the values 
of specific parameters for that fluid. These values must then 
be hard-coded. A review of the entire procedure is presented 
to aid in future lubricant additions. This review utilizes material 
from references 4 to 6. In particular, the fluid traction model 
used by the SHABERTH code is discussed in detail. 

Experimental viscosity-temperature, density, thermal 
expansion, and thermal conductivity data for RP-1 were 
obtained from reference 7. The experimental characteristics 
of traction force as a function of shear rate needed to be 
established. The traction force in a bearing contact is dependent 
not only on shear rate (or sliding velocity) but also on the 
maximum Hertzian contact pressure, the rolling velocity, and 
the temperature of the fluid in the contact. An effort was made 
to obtain experimental data on the traction coefficient as a 
function of sliding velocity for RP-1 fuel at various conditions 
of maximum contact pressure, rolling velocity, and temper- 
ature. The work was conducted on a twin disk traction tester 
by using side slip to produce the traction. Further information 
on the testing can be found in reference 8, where RP I data 
pertaining to the pressure-viscosity coefficient at ambient 
pressure as a function of temperature, which was used to 
modify the SHABERTH code, can also be found. 

Fitting the experimental data to the rheological models is 
straightforward for all but the fluid traction model. Therefore, 
only the fitting of the RP-1 fluid traction data to the fluid 
traction model will be discussed in detail. The necessary 
parameters for all of the models will be specified. 

In the process of fitting the RP- 1 traction data to the traction 
model that exists in the code, certain assumptions and simpli- 
fications were necessary. The error due to making these 
simplifications is quantified by comparing theoretical 
predictions of traction force with and without the simplifica- 
tions to each other and to the experimental data. 


Values for the parameters within the models for RP-I are listed 
in appendix A. 

Although the SHABERTH code requires the user input file 
to be in SI units and furnishes the output in SI units, the internal 
calculations are all done in English units. Also, some of the 
models used are based on empirical equations that were 
developed by using a mixture of SI and English units. 
Therefore, some variation occurs in the system of units 
required tor each piece of property data that must be hard- 
coded. The proper units for each piece of information are 
specified. All symbols are defined in appendix B. 

Density, thermal expansion coefficient , and thermal 
conductivity. The SHABERTH code assumes a linear 
relationship between fluid density and temperature, with the 
thermal expansion coefficient being the slope. The necessary 
parameters are the fluid density p (in g/cm 3 at 60 °F) and the 
thermal expansion coefficient G (in (g/em 3 )/°C), which can 
be obtained as the slope of the line of the specific gravity 
plotted as a function of temperature T. The code internally 
converts G to degrees Fahrenheit and calculates density using 

p T = p (at 60 °F) - G (at T- 60 °F) (1) 

The thermal conductivity of the fluid K/, which must be 
expressed in the code in watts/meter °C, is assumed to be 
constant with temperature. In actuality, it varies slightly with 
temperature, so the user must choose the specific temperature 
at which to take the thermal conductivity data. 

Viscosity as a function of temperature . — The kinematic 
viscosity v in centistokes (cSt) at atmospheric pressure is 
calculated as a function of temperature from Walther’s relation 
(ref. 6) 

j °gH)[log| 0 (^ + 0.6)] = A - B log 10 (T + 459.7) (2) 

A and B are lubricant-dependent constants determined by 
substituting into equation (2) the experimental data of kinematic 
viscosity v (in cSt) as a function of temperature T (in °F). 

SHABERTH calculates the absolute viscosity at ambient 
pressure r} { \ in ccntipoise (cP) using 

Vo = vp ( 3 ) 


Review and Discussion of the SHABERTH 
Rheological Models 

Fluid Property Models 

This section describes the density, thermal conductivity, 
viscosity and temperature, and pressure-viscosity and temper- 
ature models already in the SHABERTH code (see ref. 6). 


The necessary parameters for the code are A, B, and v at 
both 100 and 210 °F, in cSt. 

An auxiliary temperature-viscosity coefficient 0 (in 1/°R), 
which is needed for use in a film-thickness thermal reduction 
factor, is found by calculating 


0 = 


0.00909 In 


>)o (at 100 °F) 


(4) 


>?o (at 210 °F) 



Pressure-viscosity coefficient as a function of 
temperature. — The SHABERTH code calculates the value of 
the pressure-viscosity coefficient a (in 2 /lb) at atmospheric 
pressure and a given temperature by using the following 
relationship developed by Fresco (ref. 6): 

a = (2.303 X 10 ~ 4 ) | C F + D F logoi') + E F [log 10 (j')j 2 | " 

(5) 

where 

TF 560/(T+ 459.7) 

v v.T, 

C F ,D h ,E F constants determined by Fresco 
T temperature, °F 

The values of C F , D F , and E h are constant for all lubricant 
types, so no further data arc needed for the particular fluid. 
However, whether this model w ill adequately predict pressure- 
viscosity coefficients for fluids that differ greatly from the original 
lubricating oils in the code is questionable. This will be discussed 
in the section Pressure- Viscosity Coefficient for RP-1. 

Traction 

In this section, the traction curve is defined and the 
SHABERTH traction model is described. All symbols arc 
defined in appendix B. 

The traction curve. — Tevaarwerk (ref. 9) defines traction 
as “the ability of a fluid film, trapped under high pressure 
in the elastically deformed region of tw'o loaded curved 
elements, to transmit a tangential force from one clement to 
the other.’ ’ The characteristics of traction variation with shear 
rate for a fluid have an important effect on the behavior of 
the contacting bodies separated by the fluid, both in the motion 
of the bodies and in the amount of heat generated between 
them. When two contacting rolling elements that are separated 
by a fluid have different surface speeds, tangential forces 
develop at the area of contact. These forces, which arise from 
the shearing of the fluid layer, are a function of the rate ol 
shear of the fluid. They are also a function of the type of fluid 
and of the maximum Hertzian contact pressure, rolling 
velocity, and fluid temperature within the contact area. 
Therefore, the traction characteristics of a fluid must be 
determined experimentally for each fluid. This can be done 
by rotating two disks against one another and introducing a 
measurable amount of slip between them while using the fluid 
in question as the lubricant. Curves of traction force as a 
function of sliding rate can be obtained at various conditions 
of maximum contact pressure P ( „ rolling velocity K, and fluid 
temperature T. Figure 1 shows a rolling-element contact 
shearing a fluid film of thickness h. The variable P is the 
idealized Hertzian pressure distribution across the contact. 
Shear stresses 7 \, which develop because of the difference in 



Disk 2 


Figure l . — Rolling-element contact showing film of oil under shear; 

= £/, - U 2 * 0 . 

surface speeds, affect the heat generation and the tangential 
forces acting on the rolling elements. 

In an actual contact, the shear stress and the shear rate vary 
across the thickness of the film. It is common practice, 
however, to treat the problem as if the shear takes place at 
one plane of the Film, with the fluid entrained by disk 1 
traveling at surface speed U ] and the fluid entrained by disk 
2 traveling at surface speed U 2 . Therefore, the relationship 
of shear stress to shear rate can be characterized by measuring 
the total traction force as a function of measurable values of 
sliding speed w v , or slip, where 

u s = U x - U 2 . (6) 

The entrainment, or rolling, velocity V of the fluid is taken 
to be 

V = (7) 

2 

Figure 2 shows the typical shape of the curve for traction 
as a function of sliding speed. Traction force is expressed in 
terms of a traction coefficient p w here 


traction force 
normal load 


( 8 ) 


Three distinct regions can be identified on this curve. The 
initial low-slip region of the curve is linear and is thought to 
be isothermal in nature. At some sliding speed the traction 
behavior becomes nonlinear, though still increasing. This region 
is also thought to be isothermal. In the third region, the heat 
generated by dissipative shearing of the fluid is no longer 
negligible. This region is characterized by either decreasing 
traction with further increases in sliding velocity or a flattening 
of the curve towards a horizontal asymptotic line. 

The magnitude of a traction curve at particular conditions 
can be characterized by the maximum traction coefficient 
reached on the curve. Again, this is a function of maximum 
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f igure 2.— Traction coefficient n as a function of sliding velocity u K . 

contact pressure, rolling velocity, and fluid temperature. The 
maximum traction coefficient is denoted by p* . and the sliding 
velocity at which p* occurs is denoted by u* 

The SHABERTH traction model. — This discussion applies 
to the SHABERTH/SKF traction model, not the 
SHABFRTH/NASA model. The differences between the two 
are discussed in reference 4. The complete SHABERTH/SKF 
traction model comprises an asperity traction model and a 
fluid-film traction model. The lubricant traction characteristics 
are required only for the fluid traction model. SHABERTH 
calculates h/o. that is, the ratio of film thickness h to composite 
surface roughness o. For h/a < 0.4, traction is modeled as 
purely asperity contact. For h/o > 3.0, the traction model is 
purely lubricant-dependent. For 0.4 < hfa < 3.0, the model 
is a combination of the asperity and fluid-film models (ref 6). 

The development of the SHABERTH fluid traction model 
is explained in reference 5. Since English-system units are used 
in the reference material, they are used here to describe the 
model and to fit data to it. Basically, the energy and momentum 
equations were developed in order to find the governing 
dimensionless terms. Such terms for heat conduction and 
convection, slide-to-roll ratio, and the speed-viscosity product 
(U = r} () Va/RJ are included in these equations. The total 
thermal effect was assumed to be a multiplicative power 
function of the individual thermal dimensionless terms, for 
which the exponents were unknown. Adding a viscoelastic 
correction to the effective viscosity term produced a 
relationship between a nondimensional traction coefficient term 
p' and a nondimensional sliding-speed term (Further 
explanation of the nondimensionalization process can be found 
in ref. 5.) The exponents of the thermal terms were determined 
by using Johnson and Cameron’s experimental data for a 
mineral oil. Shell Turbo 33 (ref. 10). The ensuing curve of 
^ ' as a function of 0, 1 , 2 is constant over all maximum contact 
pressures, rolling velocities, and temperatures, for all choices 
of lubricants; it is the curve to which data for all other lubricant 
choices for the program are fit. The maximum nondimensional 


traction coefficient term on the curve is designated by (p')\ 
and (i/'d “)* is the nondimensional sliding-speed term at which 
(p') occurs. These values which were obtained by using the 
Johnson-Camcron data, arc constant: 

ip') = 0.235 and (i^j °) = 8.75 X 10 ° (9) 

From these values an equation is obtained for the maximum 
traction coefficient p* at particular conditions, and for the 
associated sliding speed u* as a function of those conditions 
(ref. 5), as follows: 


0.235=^* 


( hK f \ o: Y Kf V 

\PjCjWr) V r) v j3V 2 J 



and 


( 10 ) 


8.75X 10 


"A / hK f y 2 y^^ 2 Y 4 Aq^Y 7 
v)\ PfCjVfr) V Kf ) \ *v / 

( 11 ) 


Four other lubricants have been added to the SHABERTH 
code. They are a polyphenyl ether (5P4E), MIL-L-7808 oil, 
MIL-L-23699 oil, and a fluorinated polyether (Freon E-l ). 

Calculation of the Fluid Traction Parameters for 
SHABERTH 

In this section, the procedure for finding the fluid traction 
parameters is described. Thereafter, a discussion of concerns 
with the model and with the potential difficulties that can be 
encountered in fitting the data is presented. 

There are three basic steps involved in fitting fluid traction 
data to the SHABERTH fluid traction model: (1) traction data 
obtained for an elliptical contact must first be transformed into 
equivalent traction data that would occur over a rectangular 
contact at the same conditions; (2) the parameters that govern 
the values of p* and u* must be determined from the trans- 
formed values of the experimental data; and (3) parameters 
governing the shape of the curve of p as a function of u s are 
determined from the transformed data. 

Transformation from elliptical - to line-contact data.— The 
SHABERTH traction model was developed on the basis of a 
rectangular, or ‘dine,” contact (ref. 5). When analyzing an 
elliptical contact, the SHABERTH code divides the contact 
into discrete rectangular strips, the length of the strips being 
in the rolling (x) direction as in figure 3. SHABERTH then 
calculates a traction force for each strip and integrates the 
traction force over the elliptical contact. Experimental traction 
data obtained over an elliptical contact must be transformed 
to the equivalent traction data that would occur over a line 
contact at the same contact conditions. McCool et al. (ref. 5) 
developed a method for this transformation. 
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Figure 3 . — Rolling-element contact ellipse showing discrete rectangular strips 
utilized by SHABKRTH. 

The pressure distribution over a line contact is assumed to 
be constant in the rolling direction and to vary semielliptically 
(as shown in fig. 1 ) along the semimajor ( v) axis; that is, strip 
to strip. If the pressure distribution and geometry of the 
equivalent line contact, and the traction force-pressure 
relationship are known, the traction force over the contact can 
be calculated. The traction force-pressure relationship, which 
varies with rolling velocity, temperature, and sliding velocity, 
can be found from the experimental elliptical traction data. 

The measured total traction force T R is divided by the 
length 2 a of the elliptical contact in the y-direction to obtain 
an average traction force per unit length T R . Then a 
relationship must be found between T R and P () \ this is possible 
with the assumption that a polynomial relationship can 
reasonably model the data for T R as a function of P ( „ while 
remaining easy to integrate analytically. Therefore, a general 
relationship ot the form 

j H = d () + d\P 0 + d 2 P; } + dyPl + d 4 P A () ( 12 ) 

where 

r/ 0 ,r/|,</i, . • • polynomial coefficients to be determined 
T r T r !2o 

j r measured traction force 

a contact semimajor length 

was integrated over the line contact to obtain an expression 
for equivalent average traction per unit length, over a line 
contact, as a function of P 0 . From this integrated expression 


and the original polynomial of equation ( 12 ), the following 
relationship between the elliptical-and the line-contact traction 
coefficients was obtained (it requires only the polynomial 
coefficients and the maximum Hertzian pressure 

of the particular elliptical contact): 

_ { 4 ) G 0 + G l P„ + G 2 P; + G i P’ + G4 rt ( , 3) 
Mclliptical V3 tT/ d {) + d\P {) T d 2 P~, "1" dyP 0 T d 4 P () 

where 

G 0 (7t/2)^ 0 
G) 2d , 

G 2 (37t/4 )d 2 
G } (8/3 )<h 
G 4 (15tt/16 )d 4 

For a particular fluid, at each condition of rolling velocity, 
temperature, and sliding velocity, the experimental data for 
average traction force T R as a function of maximum contact 
pressure P (> are fit to equation (12) in order to find the 
polynomial coefficients. These coefficients and the particular 
contact pressure of interest are substituted into equation (13), 
whereby equivalent traction coefficient data over a line contact 
are determined. 

Traction parameters governing the maximum traction 
coefficient and the sliding velocity of maximum traction.- 
The traction parameters needed by SHABERTH to find the 
values of p* and u * at particular contact conditions are those 
dealing with the viscoelastic relationship and some 
proportionality constants. 

The viscoelastic relationship defines an effective fluid 
viscosity rj c , at contact conditions, as a function of ambient 
fluid viscosity r? () , contact pressure P () , and fluid entrainment 
velocity V. The viscoelastic model in SHABERTH is of the form 



where / P a and (K/K 0 ) -X{) are viscoelastic corrections for the 
pressure and the entrainment velocity, respectively. The 
viscoelastic parameters that must be found for each new 
lubricant are / P„ and X<). 

The exponent \) is obtained by fitting the experimental data 
for p* as a function of K, at constant temperature and pressure, 
to a relationship of the form 

f = c m v y > ( |5) 

where C„, is a constant of proportionality and by defining 7 , as 

7 , = « 3 + 0.5Xo(«, + a 2 ) ~ \) ~ <’1 + ”i («2 - 0 + °-5«2 

( 16 ) 



where 


— 0.236 
a 2 = 0.55 
a } = 0.22 
<*, = 0.022 


The coefficients were developed from the 

Johnson-Cameron data. The term /*, is the exponent of the 
rolling speed from the film thickness relationship. With a value 
of n j =0.7 taken for a typical unstarved contact, the equation 
reduces to 


0.165 - 73 
0.607 


(17) 


The function/ P (} is found by using equation (10) and the 
experimental contact conditions and traction data to calculate 
values of r\ e as a function of P () . From the viscoelastic 
relationship, 


f\Po = 



(14) 


f P () can be found as a function of P (t since Xq has been 
calculated previously. The constant V 0 will be discussed shortly. 

The SHABERTH model for f\P 0 ] was developed from the 
Johnson-Cameron data. Figure 4 shows the relationship 
obtained from these data plotted on log-log coordinates. This 
figure shows not only that the relationship is exponential but 
also that the value of the exponent changes at a certain 
pressure. Therefore, the relationship is modeled as a function 
of the form 



where 

A, = At , for P tt < P, 

A t = A 2 * for P () > P, 

and .4,, P,, and C„ are lubricant-dependent constants. 


2x10 6 



Figure 4. -Variation of viscoelastic function f\PJ with maximum contact 
pressure for Shell Turbo 33 oil used to develop SHABERTH Hu id traction 
model. 


relating /x* and u* to P m V, ^ 0 , and /?, are obtained: 


=(C 1 )/ > “ 1 - 14 [-2 


0 . 61 / 1 : 


Vo 


0 . 59 ^ 0 . 48 - 0 . 61 ^- 0.45 (J9) 


and 


- 0 . 4 / 1 / 


= (C 2 )P<7 i — V° 4X ° - 0 . 09 ^- 0.55 (20) 

\ p \/ 


where 


The value of P , is chosen as the pressure (in lb/in? X 10 ~ 5 ) 
at which the slope of the log-log plot of f\P 0 \ as a function 
of P changes. The value of A x is chosen from the best-fit 
value of the exponent A t over all conditions prior to P,, and 
A i is the best-fit exponent at pressures greater than or equal 
to P|. The value of C n need not be determined from this data 
fit because it is combined w ith other constants to form general 
proportionality constants C, and C 2 , which are determined 
after all of the specific parameters are found. 

By substituting the viscoelastic relationship into equations 
(10) and (11) and by grouping the fluid property constants into 
general proportionality constants, the following equations, 


A t —A x , for P n < Pj 
Ai = A 2y for P (t > P, 

By substituting experimental values of fi* and u*, at the 
specific test conditions, into equations (19) and (20), best-fit 
values of C\ and C 2 can be determined. These equations are 
used in the SHABERTH code, along with lubricant-specific 
values of Xq, P h A u A 2 , Cj, and C 2 , to predict values of /x* 
and u*. 

The constant V 0 , which serves as a rolling velocity normal- 
ization factor, is not a required parameter in the SHABERTH 
code, ft is used only when a relationship between f \P 0 and 
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P 0 is fitted by using equation (14). For each of the oils dis- 
cussed in reference 5, the value chosen for K () was that of the 
lowest rolling velocity condition in the particular test matrix. 
However, for the model of f)Pf shown in equation (18) (an 
exponential relationship), the value of Vq has no effect on the 
relevant parameters that are being determined (i.e., the values 
of A u A 2 , and P,); its only effect is to shift the curve of 
log [f P () \ as a function of \og[P ( ,\ by a constant factor. The 
values of the proportionality constants C ] and C 2 are 
determined from equations (19) and (20), which do not include 
the constant K 0 . 

Modeling the shape of the traction curve .— Once the 
theoretical value of p* has been determined by SHABERTH 
for certain contact conditions, the shape of the curve of the trac- 
tion coefficient as a function of sliding velocity is needed to 
determine p at a particular sliding velocity. Experimentally 
derived traction curves are of two basic types (see fig. 5). For 
both types, p increases linearly at low sliding speeds and then 
becomes nonlinear, increasing at a slower rate. The curve shown 
in figure 5(a) reaches the value of p * and then slowly decreases 
with increasing sliding speed. The curve shown in figure 5(b) 
increases continuously to an asymptotic value of p* . SHABERTH 
uses the continuously increasing, asymptotic model; the other 
model created convergence problems in the program. 

Values must be found for the coordinates where the linearly 
increasing portion of the asymptotic curve ends. These values 
are dependent on the contact conditions. However, the 
experimental data for all contact conditions falls on a single 
curve of (p/p*) as a function of (uj u * ) ■ Therefore, 
SHABERTH models the shape of the traction curve on the 



Figure 6.— Nondimensionali/ed traction curve showing coordinates (A',,. Y fi ) 
used to define shape of curve. 


basis of one set of coordinates, (X R ,Y R ), where the lineai 
portion of the experimental curve of (p/p*) as a function ot 
(uju*) ends (see fig. 6). The SHABERTH code then 
calculates values of p using 



0 < — ^ < X B 

it* 


and 





YpX l} \ 
Yi>Xb + YfiXp ) 


— >X„ (21) 

It* 




(a) Curve exhibiting peak behavior. 

<b) Asymptotic behavior model used in SHABF.RTH code. 

Figure 5 .— Two types of behavior exhibited by curves of experimental traction 
coefficient as a function of sliding velocity. 


where 

Y n 1 - Yh 
X„ (uju*) ~ X„ 

Issues Concerning the SHABERTH Fluid 1 raction Model 

Issues relevant to the SHABERTH fluid traction model and 
the problems that can arise during the fitting ot the 
experimental traction data and the extraction of the parameters 
are discussed. 

Transformation from elliptical- to line-contact data.—C are 
must be taken when fitting a polynomial to the elliptical-contact 
data of T r as a function of P (r As previously discussed, the 
resulting polynomial models the relationship between T R and 
P as it varies over a rectangular (line) contact. Equation ( 13) 
was developed by integrating a genera! form ot the polynomial 
across incremental strips of a line-contact area tor which the 
pressure is P (l at the center strip and decreases semielliptically 
toward the edge of the contact in the v-direction. therefore 
the polynomial for each set ot contact conditions must be valid 
down to pressures below P r 

Figure 7 shows a possible example of a polynomial fit to 
experimental data. The points A, B, C, and D represent 
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Pressure, P Q , psi 

Figure 7. — hxumple ol ,i polynomial lit through four points representing 
experimental data. 


experimental data points. As P decreases, the polynomial 
reaches a minimum and then begins increasing at values 
of P not far below the value of P () for conditions at point A— 
behavior which is not valid for traction as a function of 
pressure. 

If the polynomial coefficients obtained from the original four 
data points are used in equation (13) for the maximum pressure 
condition at point A, and to a lesser degree at conditions of 
point B . the resulting Mimc^ciiipncai values are often 
implausible. This problem was first encountered with the RP-1 
experimental traction data. Since there were only four 
conditions of pressure for the fit, second-order polynomials 
were used. Reference 5 contains experimental elliptical- 
contact data and equivalent transformed line-contact data for 
a polyphcnyl ether (5P4E), which had a test matrix containing 
three conditions of pressure. When the transformation method 
was applied to this elliptical data, also using second-order 
polynomials, the same problems occurred as with the RP-1 
data. The transformed line-contact data for 5P4E, obtained 
at the lowest two conditions of P (f by using polynomials fit 
to the listed experimental traction data, often did not match 
the original transformed data listed in reference 5. Also, they 
were not reasonably expected values. When the polynomials 
for these cases were examined, they exhibited behavior such 
as shown in figure 7. The erroneous increase of T R at 
pressures still well within the contact ellipse has a profound 
effect on the resulting Miine^diipticai values. 

A method was found to solve this problem. First, fit a 
geometric relationship to the data of f R as a function of P () 
such that 


Tr = b^") (22) 


where b, and b 2 are coefficients of the geometric fit. 

An advantage of this relationship is that T R goes to zero as 
P () goes to zero, which is in keeping with physical reality. By 
using equation (22), values of T R at pressures below the 
lowest test matrix P (i were predicted. Then the original and 
piedicted data for T R as a function of P 0 were used together 
to fit a polynomial of the form of equation (12). The resulting 
set of polynomial coefficients (rf 0 , etc.) were used in 
equation (13) to obtain final values of Mtin^Meiiipncai- For 5P4E 
the results of this method were in good agreement with the 
original transformed line-contact data listed in reference 5 for 
all conditions. Another method was also tried, using an 
exponential relationship of the form T R = Z> 3 (p,f 4 ) in place 
of equation (22), but these results did not agree as well with 
the original 5P4E data. Therefore, the method that uses 
equation (22) is recommended. 

Viscoelastic parameters.— The model of/ P 0 as a function 
of P () shown in equation (18) was developed from 
experimental data that exhibited a change in the value of A t 
at a pressure P, within the experimental test matrix; it is 
embedded in the final traction equations used in the 
SHABERTH code. Therefore, values for A h A 2 , and P , must 
be obtained for each new lubricant. 

However, it is not clear what procedure should be used to 
determine values for P, and A 2 when they do not occur within 
the pressure range over which the experimental data were 
taken. The / P () data shown in reference 5 for the 5P4E and 
the MIL-L-7808 oils indicate that in neither case is there a 
change in the value of A t at any place in the experimental 
pi ess u re range. How ever, values of P\ and A 2 were selected. 

On the basis of the behavior of other fluids under high 
pressure (refs. 10 and 1 1), it is reasonable to assume that there 
is a pressure P, at which the value of the exponent A, 
decreases from A } to some A 2 . The data in reference 1 1 show' 
that the pressure-viscosity coefficient a, normally associated 
with the simple power law v = r} 0 e n ' p <\ decreases to 
P' < < a at some pressure. The pressure-viscosity model of 
equation (18) is somewhat different than this one, but since 
equation (18) is still a power law, it would also be expected 
to exhibit a decrease in the value of the exponent A, at some 
pressure. 

It is impossible to know the exact value of Pj without 
experimental data. A value should be chosen that gives some 
weight to the experimentally derived value of A x \ that is, a 
value which extends, somewhat, the range at which A\ is 
valid beyond the highest pressure in the test matrix. However, 
the value chosen for P } should not be too much higher than 
the experimental range, since the error due to using A l at 
pressures where the actual behavior follows A 2 has 
exponential effects on the predicted values of viscosity. For 
the 5F4E and the MIL-L-7808 oils respectively, values of 
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Pj were chosen that were 14 and 20 percent greater than the 
highest experimental contact pressure. This seems reasonable. 

The only data available on which to base a choice ol the 
value of A~> are those ol Johnson and Cameron for Shell 
Turbo 33. The most reasonable recourse seems to be to choose 
a value of A->IA\ based on A 2 IA\ for Shell Turbo 33 and then 
find A, from the value of A ] for the particular fluid. For Shell 
Turbo 33, A 2 IA\ = 0.635. 

In choosing a value ol A 2 !A\ lor other fluids, a more 
conservative prediction of viscosity is obtained by erring on the 
low side of the actual A 2 than by erring by the same amount 
on the high side, because ol the exponential effect. Therefore, 
a value of less than 0.635 is recommended for A 2 IA V 

For the other original lubricants in the SHABERTH code, 
the following A 2 IA X values were used: 5P4E, 0.456; 
MIL-L-7808, 0.363; and M1L-L-23699, 0.5. 

Final fluid traction equations ,— Equations (19) and (20) 
show the final equations used by the SHABERTH code to 
predict g* and u* as functions of P ( „ 170, V , and h. However, 
a slight error has been detected in these formulas. In the 
original traction equations developed by McCool, et al. (ref. 5) 


and 


M.v 


= (c 2 )P t ; 



-0.4/1, 

y h) 1 1 f <( )'4 () 09 Vr 05 ' 


(27) 


where 


Aj = A u for P (> < Pi 
Aj = A 2 , for P tt > P| 


The exponent of P (> cannot be easily changed in the 
program for the following reason: once A | , A 2 , X 0 , and P\ 
have been determined from equations (17) and (18), the 
traction data are used to determine the values of Cj and C 2 
for each lubricant. The hard-coded values of Cj and C 2 for 
the existing choices were calculated by using Pf x 14 . Without 
the original traction data for these lubricants, the values of 
C, and C 2 associated with P" 1 ' 28 cannot be calculated. The 
difference in traction prediction between the two sets of 
equations for the RP-1 data at the experimental conditions will 
be shown to be negligible. 


^ = c m P' 1 /7- (,28 /r°- 45 (23) 

where 

C m proportionality constant 
b 27t/? a (P„)/E' 

h C t/ (a^ 0 F)° 7 P r r° 3 for typical unstarved contact 


Development of the SHABERTH Fluid 
Traction Parameters for RP-1 

In this section, the experimental traction data for RP-1, a 
modified model of the pressure-viscosity coefficient for RP-1, 
and the fitting of the data to the SHABERTH fluid traction 
model arc discussed. 


Substituting the expression tor b into equation (23) gives the 
relation 

= C m P„ * ■-**/* — ,,-4S (24) 

The - 1 .28 value of the P tt exponent is that value for which 
h is still a separate variable in the equation. To find the total 
relationship between /x* and P () , the film thickness term h was 
broken up to give 


Here, the -1.14 value of the P 0 exponent is that for which 
h is no longer a variable in the equation. 

In equations (19) and (20), however, h has been reintroduced 
as a variable. The exponents associated with q {) and V are the 
appropriate ones since the film thickness term includes 
contributions of 17 q and V. The exponent of P on the other 
hand, has been left as —1.14. 

The corrected equations should read 

^ ' 1 ',n.»K <0 - 48 - aw ^ ' a4S (26) 


Experimental Traction Data for RP-1 

After a description of the RP-1 traction test matrix, the effect 
of transforming the experimental elliptical-contact traction data 
to equivalent line-contact traction data is shown. General trends 
of the transformed RP-1 data arc also presented. 

Test matrix conditions,— A twin disk tester generated data 
for traction force as a function of sliding speed for RP-1 fuel. 
An elliptical contact under various conditions of Hertzian 
pressure, rolling speed, and temperature (ref. 8) was used. 
The transverse radius of curvature of the lower disk was 
infinity; the upper disk was a toroid, the curvature of which 
produced the desired contact geometry. The traction data were 
generated with side slip, which can be produced by skewing 
the toroid about the normal to the horizontal plane. The slide- 
to-roll ratio uJV can be controlled and measured via the skew 
angle. The test matrix is listed in table I. Figure 8 shows a 
typical set of curves for experimental traction coefficient ft 
as a function of uJV tor elliptical contact. This particular set 
of curves shows the variation of traction with rolling velocity 
at one condition of pressure and temperature. In general, 
traction decreases with increasing rolling velocity, increases 
with increasing pressure, and decreases with increasing 
temperature. 
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TABli: I. TARGHT THST MATRIX CONDITION'S" 


Herl/ian pressure, GPa < kpsi ) 1.01 ( 147) 

1.27 (185) 
1.60 (233) 
1.92 (279) 


Rolling velocity. m sec (in., sec) 10 (390) 

30 (1200) 
50 (2000) 

Temperature. T i°F) 40 (104) 

65 (149 ) 

‘‘All parameters uere texted. 



Sliding velocity/rolling velocity, u s /V 

Figure 8.-F,\perimental traction coefficient as a function of slide-to-rol! ratio 
lor RP- I at varying values of rolling velocity. Contact pressure, 279 kpsi; 
temperature. 1 1 1 °F. 

The conditions shown in table I represent the maximum 
capabilities ot the test rig. The upper limits extend above 
those of the original lubricants. For example, the Johnson- 
Cameron data had the highest Hertzian pressure condition at 
225 000 psi (F nm = 1000 in. /sec), and the MIL-L-7808 
data included the highest rolling speed at 1820 in /sec 
(P„ max 150 000 psi). Values of Hertzian stress greater 
than the maximum in the RP-1 matrix would result in h !o 
values that would trigger the asperity-contact-only traction 
model. Therefore, the test matrix shown in table I is quite 
adequate for the purpose ot modeling the traction behavior 
of RP-1. 

Transformation from elliptical-to line-contact data.— The 
RP-1 data were transformed from ellipticai-to line-contact data 
by the method previously described. Since there were only 
four conditions of pressure in the fit, second-order polynomials 
were used. Figure 9 shows a comparison of the RP-1 curve 
of p as a function ot uJV tor elliptical contact to that for the 
equivalent transformed line contact for RP-1 at one set of 
conditions. The elliptical-and line-contact curves at other 
conditions compare similarly. The equivalent line-contact data 
always exceed the elliptical-contact data, the greatest amount 
of difference being in the nonlinear region of low uJV. The 
linc-to-elliptical ratio then decreases at increasing uj V. Also, 
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Iigure 9.— Comparison of RP- 1 experimental traction data for elliptical contact 
to analytically transformed traction data for equivalent rectangular contact. 
Contact pressure. 233 kpsi: velocity. 12(H) in. /sec: temperature. I(M T 


the value of u * , the sliding velocity at which p* occurs, tends 
to be lower for the line-contact data than for the elliptical data. 
The behavior ot the transformed data is typical of that for 
other oils. 

General trends of the RP-1 data.— Figure 10 shows a plot 
of p/p* as a function of uju* for RP-1 for all conditions in 
the test matrix. Since the SHABERTH model of the traction 
curve is based on the assumption that p/p * is a function of 
l, J u s °nly, it is desirable for the experimental data for p/p* 
as a lunction of u s /u* to fall on one curve. In figure 10 the 
RP-1 data all follow the same general shape, though some 
scatter is present. The scatter is not deemed large enough to 
prohibit the use of the SHABERTH traction model. 

Figure 1 1 shows the effect of maximum contact pressure 
on the values of p and u* for RP-1 at one condition of 
rolling velocity and fluid temperature. The behavior shown 
is typical over the range of conditions. The value of p* 
increases with increasing contact pressure, and the value of 
u* decreases. In figure 12. the effect of rolling velocity on 
the values of p* and u* for RP-1 is shown at one condition 
of contact pressure and temperature. Again, the behavior 
shown is typical over the range of conditions. The value of 
P* decreases with increasing rolling velocity, and the value 
of u* increases. 

Pressure-Viscosity Coefficient for RP-1 

The SHABERTH code calculates a pressure-viscosity value 
as a function of temperature and lubricant ambient kinematic- 
viscosity by using an expression developed by Fresco (ref. 6). 
This semicmpirical model was developed by using data for 
fluids other than RP-1 , and it accounts for different lubricant 
types only in the value of the viscosity. The coefficients used 
in the calculation are constant for all lubricant types. This 
model was not used in the analysis of the RP-1 data. 
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Nondi memorialized sliding velocity, u s /u s 

. . , , , function of u /«* for all test matrix conditions of maximum contact pressure, rolling velocity, and 

Figure 10. — RP- 1 experimental data tor lit fi as a Junction oi u x /u x ior an 

fluid temperature. 



140 1 60 1 80 200 220 240 260 280 

Maximum contact pressure, P 0 , kpsi 

ijgure 1 1 — Variation of RP-1 experimental maximum traction coefficient 
and associated sliding velocity with maximum contact pressure. Velocity, 
2000 in. /sec; temperature, 104 °F. 



400 800 1200 IbUU 

Rolling velocity, V, inJsec 

Figure 12. Variation of RP-1 experimental maximum traction coefficient 

and associated sliding velocity with rolling velocity. Pressure. 233 kpsi; 
temperature, 104 D F. 


As p;ir( of the effort to measure and analyze the traction 
data for RP- 1 . pressure-viscosity data generated by Bridgman 
for kerosene were obtained (ref. 12). The traction data specific 
to RP- 1 were also used to check the validity of Bridgman's 
pressure-viscosity data for use with RP- 1 . By using the Barus 
prcssiire-viscosiiv relationship, that is. 


where 


}) — y 0 exp 



(28) 


v] 6 ambient absolute viscosity at inlet temperature. Pa-see 
0 inlet temperature, °C 

D r lubricant-dependent pressure solidification 
temperature, °C 

B r lubricant-dependent constant. °C/Pa 
P pressure. Pa 

and the reference 8 experimental data, the constants B r and D, 
were determined for RP-1 as B r = 544 °C/Pa and D p = 25.6 °C 
(rel. 8). This was the pressure-viscosity calculation used in 
fitting the RP-I traction data to the SHABERTH traction model. 

The pressure-viscosity coefficient a. which is used by 
SHABI.R I H, can be calculated from the Barus model bv usins 

d + D r K ’ 

and converting from units of Pascals 1 to square inches per 
pound. 

Figure 13 shows the RP-I pressure-viscosity coefficient as 
a function of temperature for both the Fresco and Barus 
models. The Fresco model predicts a nearly linear relationship 
between the pressure-viscosity coefficient and temperature, 
whereas the Barus model predicts a hyperbolic relationship. 
The difference between the two models' pressure-viscosity 
values w<as deemed significant enough to justify changing the 
SHABERTH code so that it would branch to equation (29) 
(with the appropriate values of B, and D p ) instead of using 
the Fresco model when the lubricant is RP-1. 


Fluid Traction Parameters 

The final values of the SHABERTH fluid traction parameters 
obtained by fitting the RP-I traction data to the previously 
discussed models are listed in appendix A. Also specified arc 
other changes that need to be made to the SHABERTH code 
to run it with RP-1 as the lubricant choice. 

Viscoelastic parameters. -The value of X () was determined 
from equation (17) alter a best-fit value of y, was found. 
This value was obtained by varying y, from its lowest to 
highest experimentally derived value, while comparing the 
predicted p to the transformed experimental p over all 
conditions, and choosing the y, that gave the lowest 



figure It Variation of pressure-viscosity coefficient for RP-I with 
temperature— comparison of I- reset, model (lluid-independcnt) and Barus- 
Vogel model. 


average percent error over the pressure and temperature 
range. The resulting value was y, final = -0.28. which gives 
To final = 0.733. 

The function/ P 0 was calculated from equations (10) and 
(14) by using X (l = 0.733 and V 0 = 390 in. /sec. The 
relationship between f P 0 and P„ is shown in figure 14 at a 
particular rolling velocity and temperature; this behavior is 
typical at all conditions. No change in slope occurs in the plot 
of lo Pl/ p » I with respect to log|/»„], thereby indicating that 
the expected decrease in the exponent A, does not occur 
within the experimental range of maximum contact pressure. 
This is true at all conditions of speed and temperature in the 
matrix. However, based on the behavior of other fluids under 
high pressures (refs. 10 and 1 1 ), an assumption can be made 
that for RP-1, there is a /*, at which the value of the exponent 
A, decreases from A, to A 2 . Since P, cannot be determined 
from the present experimental data, the value of A h which 
was obtained from experimental data, will be assumed to be 
valid up to a pressure 20 percent higher than the highest 
maximum contact pressure in the test matrix. This choice gives 
some weight to the experimentally derived data, while 
providing a safety valve against predicting unreasonably large 
values for viscosity. Therefore, /*, = 340 000 psi. 

In order to model/;/*,,; as a function of P tl only, the data 
for f P„ with respect to /»„ for all conditions of rolling 
velocity and fluid temperature should fall on the same curve 
Figure 15 shows the variation of /;/*,, with P 0 for all of the 
conditions and shows the best-fit line. Note that a considerable 
amount of scatter exists. 

Figure 16 again shows the variation of f ,P„ with P„ for 
all conditions, but the data sets taken at the two fluid 
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Figure 14. —Variation oi viscoelastic function with maximum contact pressure 
for RP 1. Velocity. 1970 in. /sec; temperature. 110 °F. 
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Figure 15. — Variation of viscoelastic' function with maximum contact pressure 
for RP-1 at all test matrix conditions. 
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Figure 16. — Data and best-fit lines showing variation of v iscoelastic function 

with maximum contact pressure tor RP- 1 at two temperatures. 

temperatures are plotted separately and the best-lit line is 
drawn through each. From this figure, it is apparent that 
temperature has an effect on the relationship. The value ol 
A x remained constant at both temperatures; only C„ is a 
function of temperature. Therefore, A x can be calculated 
without regard to the effect of temperature. The best-fit value 
obtained for A\ was 3.23. 

The value of A 2 cannot be determined directly from the 
experimental data since it occurs outside of the test matrix. 
The best alternative is to choose a value for A 2 !A X based on 
the behavior of other fluids and determine A 2 from the value 
of A i for RP-1. The only fluid data in the SHABERTH 
reference material showing a decrease in A x within the 
experimental test matrix are those for Shell Turbo 33, 
obtained by Johnson and Cameron (ref. 10). Their data give 
I ) Shell Turbo 33 = 0.635. 

For the same reasons stated in the section Issues Concerning 
the SHABERTH Fluid Traction Model, A 2 IA { for RP-1 was 
chosen as a lower value than that for Shell Turbo 33. For 
convenience, the value chosen was (/E/^iW-i = ^-5, giving 
A 2 = 1.6. 

An additional point can be made here concerning the RP- 1 
parameters. The particular value chosen for A 2 is not of 
critical importance to the traction prediction for RP-1. The 
range of pressures in which A\ is valid should cover all cases 
in which the fluid traction plays a significant role. At higher 
contact pressures, the traction force prediction is governed 
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TABLE II.— VALUES OF C, AND C 2 FOR VARIOUS CASES 


Equations 

used 

Tom pc rat u re - a ve raged 
values 

Temperature-dependent values 

At 104 °F 

At 149 °F 

Cl 

C, 

Cl 

C 2 

Cl 

c 2 

(19) and (20) 
(26) and (27) 

1 13 450 
629 519 

24.65 
1 36.40 

126 456 
700 536 

2 1 .40 
118.01 

99 262 
552 045 

28.19 

156.46 


predominantly by the dry asperity traction model. However, 
a value tor A 2 should be coded in case it is required. 

There is currently no mechanism in the SHABERTH 
program to handle a temperature-dependent f\PJ. Because 
there were only two conditions of temperature in the RP-1 
experimental matrix, this study describes the error involved 
in ignoring the temperature dependency, but it does not 
develop a temperature-dependent algorithm. 

Proportionality constants . — The temperature -dependent 
constant C n in equation (18) is contained in the general 
proportionality constants C\ and C 2 , which must be found by 
fitting the transformed experimental RP-1 p* and u* data to 
equations (19) and (20), once the viscoelastic parameters have 
been determined. Therefore, C| and C 2 reflect the tempera- 
ture dependency encountered in the f\P () data. The values 
chosen for these constants are the average of the best-fit values 
obtained at each temperature. The error resulting from the 
temperature-averaged values of C] and C 2 will be quantified 
in the section Results and Comparisons by comparing the 
traction predicted by using them to the traction predicted by 
using the values of C t and C 2 obtained at each temperature. 

The values of C ] and C 2 are also affected by the potential 
correction to equations (19) and (20), which led to equations 
(26) and (27). Therefore, in the Results and Comparisons 
section, the correction to the exponent of P () wall be shown 
to have a negligible effect on the resulting value of p* . 
Equations (19) and (20) are used in all subsequent calculations 
to determine values of C, and C 2 . The values ultimately 
chosen for use in the SHABERTH code were the temperature- 
averaged values found by using the uncorrected equations (19) 
and (20); that is, C, - 113 450 and C 2 = 24.65. 

Table II lists the values of C, and C 2 found for the various 
cases described above, that is, the temperature-averaged and 
temperature-dependent values needed for the uncorrected 
equations (19) and (20) and the corrected equations (26) and (27). 

Traction curve shape parameters . — As stated previously, 
the traction curve predicted by SHABERTH follows a fixed 
mathematical function, increasing linearly at low' sliding speeds 
and then asymptotically approaching the value of p* . Values 
of X B and Y B obtained from the nondimensionalized curve of 
p/p* as a function of u s /u* are used to calculate p as a 
function of u s from equation (21), once p* and w v * are 
known. The values of X B and Y B that w ? ere found from the 
RP-1 data were 0.15 and 0.65 respectively. 


Results and Comparisons 

In this section the effect of the correction to the exponent 
of P 0 in equations (19) and (20) will be shown to be 
negligible. The term “experimental data“ will be used to mean 
the equivalent line-contact data that was transformed from the 
experimental elliptical data. Figure 17 shows a comparison 
of curves of p as a function of u s /V. These curves were 
generated by the SHABERTH model at test matrix conditions 
by using the corrected ((19) and (20)) and uncorrected ((26) 
and (27)) equations with the appropriate values of C, and C 2 
for each. For reference, the experimental curve at the same 
conditions is also shown. For both SHABERTH curves, the 
effect of temperature on f P a has been neglected (i.e., C, and 
C 2 are temperature-averaged values). The value of the 
maximum traction coefficient of the corrected curve exceeds 
that of the uncorrected curve by 5 percent. Over the range 
of conditions, the maximum difference in the two predicted 
values of p* was 5 percent, with an average difference of 
3 percent. The percentages for maximum and average 
differences in the values of a* were comparable to those 
of p*. 

In figure 18, which show s corrected and uncorrected curves 
at an extrapolated condition of pressure (330 kpsi) and rolling 
velocity (3937 in. /sec), the curves differ by 6 percent. 

As w ? as stated previously, the values of the proportionality 
constants C ] and C 2 that were hard-coded for the original 
lubricants would need to be recalculated by using the original 
traction data if equations (26) and (27) were to be substituted 
for equations (19) and (20). However, for a I to 5-percent 
difference in the value of p*, recalculation does not seem worth 
the effort. 

Neglecting the temperature dependency of C, and C 2 , 
though, has a greater effect on the results. Figures 19 and 20 
compare experimental traction curves at various conditions, 
not only to those predicted by using temperature-averaged 
values of C\ and C 2 but also to those predicted by using 
temperature-dependent values of C { and C 2 . 

The curves in figure 19 are all at the lower temperature 
condition (104 °F). In all cases the temperature-dependent 
curves predict the value of p* quite well. The errors between 
the expet imental values of p* and those predicted by using 
temperature-dependent C, and C 2 are, respectively, 2.6, 0.6, 
0.3, and 1.2 percent for the four cases show n. Over the range 
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Traction coefficient, fj. 



Sliding velocity/rolting velocity, u s /V 


Rcine 1 7. -Comparison of SHABERTH predictions for RR- 1 using uncorrectod and corrected 
I 181 m sec: temperature, 95 °F. (Transformed experimental data shown lor reference.) 


traction equations. Contact pressure. 147 kpsi; velocity 


Traction coeffiaei 



Traction coefficient, fi 





(a) Low maximum contact pressure: P 0 , 147 kpsi (velocity. 1953 in. /sec). 

(b) High maximum contact pressure: P ir 279 kpsi (velocity, 1969 in.; sec). 

(c) Low rolling velocity: V , 445 in. /sec (P ir 185 kpsi). 

(d) High rolling velocity: K 1969 in. /sec (P n . 185 kpsi). 


Figure 19.— Traction coefficient of RP-1 as a function of slide-to-roll ratio at 104 ^-comparison of curves from transformed experimental 
from SHABERTH predictions using temperature-averaged and temperature-dependent coefficients. 


data with those 


4x10" 2 



4x10~ 2 



<ai High maximum contact pressure, high rolling velocity: P tt , 279 kpsi; V. 1953 in., see 
(hi Low maximum contact pressure: P tf . 185 kpsi (velocity, 1969 in /see) . 

(cl Low rolling velocity: V, 1200 in 'sec (P,, 279 kpsi). 


t 'f Ua ' c " c ' llticm of RP “ 1 as a func,i »» »' slidc-lo-roll ratio at 149 “F comparison ol curves Iron, transformed 

1 10,11 ^HABhRIH preJiciions using temperature-averaged and temperature -dependent coefficients. 


experimental data with those 


of conditions at the lower temperature, the average error was 
3.4 percent and the highest error, 8.2 percent. 

Since the value of p* decreases with increases in tempera- 
ture. the use of the temperature-averaged values of C, and C\ 
underpredicts the experimental values at the lower 
temperature. In the cases shown, the values of>* from using 
temperature-averaged C \ and C\ differ from the experimental 
ones by 8.2, 9.9, 10.0, and 9.3 percent respectively. The average 
error in fx* over the range of conditions at the lower tempera- 
ture was 10.3 percent, and the highest error was 17.9 percent. 

Figure 19 also shows that the accuracy with which the shape 
of the traction curve is predicted varies with contact pressure 
and rolling velocity. As a result of SHABERTH s asymptot- 
ically increasing model of the traction curve, any thermal 
dissipation effects, that is, decreasing traction at high it V /K, 
ate completely neglected. Of greater importance, however, 


is the prediction of the low uJV region, since ball sliding in 
a bearing occurs predominantly at low itjV. 

Parts (a) and (b) of figure 19 show' curves at conditions of 
different maximum contact pressures, but of the same rolling 
velocity. At low contact pressure (fig. 19(a)) the effect of 
thermal dissipation is small, and the predicted temperature- 
dependent curve agrees well with the experimental one at high 
uJV. At low ujV s however, the SHABERTH curve both 
overpredicts the initial slope and remains linear up to a higher 
uJV than does the experimental curve. At high contact 
pressure (fig. 19(b)) the thermal effect at high uj E is much 
more pronounced and, thus, not well predicted by 
SHABERTH. However, the initial portion of the curve is better 
predicted at the higher contact pressure, almost exactly 
matching the shape of the experimental curve up to 
0.01 5 uJV. 


is 



Parts (c) and (d) of figure 19 show the effect of rolling 
velocity on the shape of the curve at the same contact pressuie. 

At low rolling velocity (fig. 19(c)) the high ujV region is 
well predicted because of the absence of significant thermal 
effects, but the initial slope is underpredicted. Also, the 
predicted curve 1 villains linear up to a higher » ( /Ethan does 
the experimental curve. The high rolling velocity curves in 
figure 17(d) show that the thermal effect increases with rolling 
velocity, thereby causing greater discrepancies at high ujV. 
The initial portion of the curve is, however, very well 
predicted, accurately matching the shape of the curve up to 
0.06 ujV. 

The results at the higher temperature (149 °F) were very 
similar (see fie. 20). Again, the values of ft* obtained by using 
temperature-dependent C, and C : show good agreement with 
the experimental ones. For the conditions shown in figure 20 
parts (a) to (c), the temperature-dependent prediction of ft 
differs from the experimental prediction by 5.9, 1 .5. and 3.9 
percent respectively. The average error over the range of 
conditions at this temperature was 4. 1 percent, and the highest 
error. 9.8 percent. The values of ft* determined by using 
temperature-averaged values of C\ and C> overpredicted the 
experimental n* by 21 . 1 . 12.2. and 9.9 percent in figure 20 
parts (a) to (c). respectively. The average error was 14.4 
percent and the highest. 25.3 percent. 

The same trends in ability to predict the shape of the traction 
curve are seen at the higher temperature as at the lowei . 
Figure 20(a) shows curves at a condition of high contact 
pressure and high rolling velocity. The initial portion of the 
curve is very well predicted up to 0.03 u s /V, whereas the 
hi ah ujV region is not well predicted because of thermal 
dissipation effects. At a lower contact pressure and the same 
rolling velocity (fig 19(b)). the high u s /V region is well 
predicted, but the predicted curve at low ujV remains linear 
up to a higher ujV ratio than does the experimental one. 
Figure 20(e) shows curves at the same contact pressure as in 
figure 20(a) but at a lower rolling velocity. As in the low- 
temperature ease, the initial slope of the curve is 
underpredicted. There is still a significant thermal effect in 
this experimental curve because of the high contact pressure. 

In general, the error from using temperature-averaged values 
of C'i and C\ runs approximately 10 percent higher than that 
from using the discrete, temperature-dependent values. The 
traction model for RP-1 does very well (1) at predicting the 
shape of the low ujV region of the traction curve at 
conditions of high maximum contact pressure and high rolling 
velocity and (2) at predicting the shape of the high ujV 
recion at conditions in w hich thermal dissipation does not play 
a significant role (i.c.. low maximum contact pressures and 
low rolling velocities). At conditions of low contact pressure 
or low rolling velocity, the initial slope of the traction curve 
is not well predicted, and the predicted curves remain linear 
up to higher ujV values than do the experimental curves. 


Concluding Remarks 

An effort is underway both to provide data on and to develop 
a bearing design methodology for Rocket Propellant I (RP-D, 
a hydrocarbon fuel, for possible use in future engine programs. 
Integral to this effort is the modification of the rolling-element 
bearing analytic code SHABERTH to incorporate RP-1 as a 
bearing lubricant. Modification of the code requires that certain 
specific property data be hard-coded. Of the required data, 
the only information not previously known tor RP- 1 were the 
characteristics of traction force as a function of shear rate; 
these characteristics must be determined experimentally. The 
traction data were obtained and fit into the traction model of 
SHABERTH, thereby providing a set of parameters that were 
hard-coded. 

A slight error was detected in the traction equations that aie 
in the SHABERTH program. However, lack of traction data 
on the original lubricants prevented immediate correction of 
the error. The error was found to make an insignificant 
difference in predicting traction for RP- 1 . A particular traction 
coefficient for RP-1 was found to be temperature-dependent; 
this is not accounted for by the model. Prediction of the 
maximum traction coefficient by using temperature-averaged 
coefficients erred by an average of 12 percent, whereas 
prediction by using the discrete coefficients at the two 
temperature conditions erred by an average of 4 percent. 
Currently, the temperature-averaged coefficients are being 
used. Future work may include developing and adding a 
temperature-dependent coefficient into the model. The model 
of the pressure-viscosity coefficient as a function of 
temperature was tailored to RP-1 by using RP-1 pressure- 
viscosity data. 

A potential project is the investigation of more recent traction 
models, with the possible replacement of the present traction 
model in the SHABERTH code in mind. Much work has been 
done in the area of traction modeling since the development 
of the SHABERTH code. It is now possible to reduce the 
traction data to a few fundamental parameters from which the 
entire traction curve can be predicted with good accuracy, 
including the high slide-to-roll region where thermal 
dissipation can become a significant factor. 

One such model, which has shown very good agreement 
with experimental data on a wide variety of lubricants and 
which would be a good candidate to replace SHABERTH, is 
the Johnson-Tevaarwerk model. In addition to its ability to 
accurately predict fluid traction due to ball sliding, the Johnson- 
Tevaarwerk model has demonstrated the capability to take into 
account spin and asperity contact as well. Further experimental 
work on RP-1. as well as work on the SHABERTH code, 
would be needed for the spin and asperity options. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 17. 1990 
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Appendix A 

Changes to the SHABERTH Code 


This appendix describes the changes which must be made to the SHARPptu i • , 

d«.icc. I, includes .he modified poriious „f ,he SHABERTH c J e ' n '° ' “ “ l “ briC '"" 

The 'l!l' l |t"d‘'he,cu, > 'h‘i's all or fh~'nv“ “ ddCd f° ' h ' COde in l986 ' Somc verslons ° r the <*& ™y no. have this addilioo. 
htrun has all ot the necessary information tor adding both the Freon E-l and the RP-1 choices. 

VISCQ2 


^ £ ~ — ca,;^ Iisi in 

DD C A A / n . _ 


BB = 544. 


(constant B, in Barus pressure-viscosity relationship, l/GPa) 


DP - 25.6 (pressure solidification temperature D p in Barus pressure-viscosity relationship, °C) 

VISC02 


.SUBROUTINE VI5C02 ( NCODE, A , B, TEMP , RH06 0, G, VISCP, ALO, RHO ) 

BB = 544. 

DP = 25.6 

TT = (TEMP - 32.)x(5./9.) 


C = 10 . X*A 

EXP = C / ( TEMP + 459.7 )XXB 
VISCS = lO.xxEXP - .6 
RHO = RH060 - G x ( TEMP - 60. ) 
VISCP = RHO x VISCS 


0 CALL ALPHAO ( Bx.2, VISCS, TEMP+459.7, ALO ) 

IF (NCODE. EQ. 6) ALO = ( BB/ ( TT+DP ) )X6 . 894E-06 

12 RETURN 

13 END 


EHDSKF 




Cl = 113450. 
C2 = 24.65 
C3 = 0.0329 
C4 = 0.203 
A 1 = 3.23 
A2 = 1.6 
PI = 3.4E5 
SLB = 0.15 
FBK = 0.65 


EHDSKF 


(C,) 

(Ct) 

(0.48 -0.61 Ay) 
(0.4A„ - 0.09) 
(A|) 

(At) 

(Z 5 ,) 

(*«> 

(Yb) 


SUDROUflNE EHDSKF( ET , H, V, P0, US, N, UM ) 

DIMENSION Cl(6),C2(6),C3(6),C4(6),Al(6),A2(6),Pl(6),SLB(6),rBK(6) 
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3 DATA C1/A73A. 5, A6293., 68910., 27350., 28000., 113A50./, C2/30A 8A8, 

4 1102 58,17.31,124. 4226 ,125.012,24.65/, C3/0 . 297 , -0 . 0934, 0 . 053, 

5 2-0.0873,-0.0881, . 0329/ , C4/0 . 03 , 0 . 286 , 0 . 19 , 0 . 282, 0 .284, 0 . 203/ , Al/ 

6 33.42,4.08,3.29,6.88,6.80,3.23/, A2/2. 14,1. 48, 1.5,3. 44, 3. 40,1 6/, 

7 4pi/l . 5E5 , 1 . 7E5 > 1 . 7E5, 2 . 2E5, 2 . 0E5> 3 . 4E5/ » SLB/0 . 1 ,0 . 25, 0 . 15, 0 . 2 , 

8 50.25,0.15/, FBK/0. 65, 0.68, 0.65, 0.68, 0.68, 0.65/ 


IF( N .GT. 0 .AND. N .LT. 7 ) 


GO TO 10 


10 KT = 6 

11 100 FORMATC^O ' !!l 19 ( • *'),/' 0 AN IMPROPER LUBRICANT TYPE CODE HAS BEEN 

13 ^PASSED TO EHDSKF . EXECUTION TERMINATES '/' 0 ,119C 

14 STOP 

15 10 CONTINUE 

16 UM = 1 .E 8 

17 IF ( PO .LT. 1000. ) RETURN 

18 IF C PO .LT. P1(N) ) XI = (PO/PKN)) ** A1(N) 

IQ TFf Pfl GT P1(N) ) XI = (PO/P1CN) ) ** A2(N) 

20 UMS = C1(N> / POMl.lt * X1MO.H , ETM0.59 « VMC31N) / "«» 45 

21 USS = C2(N)/P0**0.14/Xl**0.4/ET**1.1 * V**C4(N) * H**0.55 

22 FBREAK = FBK(N) 

23 SVB = SLBCN) 

24 SV = US/USS 

25 UM = UMS * FRICTNC FBREAK, SVB, SV ) 

26 RETURN 

27 END 


I UPROP 

In the subroutine LUPROP. Text must be added to the DATA list in lines 8 and 9. In line 15, ‘,600’ must be added to the 
GOTO statement. Lines 82 to 96 must be added. These lines contain the following information. 


VIS1 = 1.571 
VIS2 = 0.7542 
A = 13.98464 
B = 5.26787 
RHO60 = 0.8065 
G = 7.2E-04 
COND = 0.13806 
BETA = 0.007278 
AKN = 20. 

FRIC = 0.070 


(kinematic viscosity at 100 F, cSt) 

(kinematic viscosity at 210 °F, cSt) 

(constant in Walther’s equation) 

(constant in Walther’s equation) 

(fluid density at 60 °F. g/cm 3 ) 

(fluid thermal expansion coefficient, 1/°C) 

(fluid thermal conductivity, W/m C) 

(auxiliary temperature-viscosity coefficient, 1/ R) 

(EHD high-contact stress factor-not used with the SKF fluid traction model) 
(Allen friction coefficient— not used with the SKF fluid traction model) 


LUPROP 

1 SUBROUTINE LUPROP ( NC0DE, KK, A, B, BETA, RH060, G, COND, VIS1 

2 $ VIS2, XLUBE, IMET ) 

3 DIMENSION AL(9,6), XLUBE(9,6) 

COMMON / FLMDAT/ EM0DC2), AKN, FRIC 


4 


21 


V • 


1 V 


5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


',' 33 ',' 

1 PL f , 1 YP ' , , HN I , 1 L V 


DATA AL /' 1 , * SH f , * EL ' , * L ' , ' TU ' , ’ RB ' , • 0 

1 1 IL 1 , T -L 1 , ' -7 • , 1 80 f , '8G* , 1 », f 1 , 1 PL 1 , 1 YP ' , 1 HN • , 1 L V'ET 1 

2 *H 1 , 1 MC 1 , 1 S2 1 , f 93 f , 1 M 1 , 1 1 L 1 , * -L * , 9 -2 * , f 36 1 , f 99 1 , 1 f J 

^ >' 1 , % SP 1 , 1 EC * , 1 1 A • , 1 L l , t El l , f *, f 1 , . , RQ* 1 CK 1 

4 'ET',' P','R0',’PE','LI 


'NT' 


1 f / 


... FILL XLUBE WITH THE APPROPRIATE ALPHMERIC LUBRICANT TYPE DATA 
DO 101=1, 9 

10 XL UBE( I , KK ) = AL ( I , NCODE) 

KT = 6 

... DETERMINE THE LUBRICANT TYPE FROM NCODE 

GO TO C 100, 200, 300, 400, 500, 600 ), NCODE 
100 CONTINUE 

. . . THE LUBRICANT IS SHELL TURBO 33 
VIS1 = 64. 

VIS2 = 8. 


20 

A = 10.34871 

21 

B = 3.672531 

22 

RH060 = 0.879975 

23 

G = 6 . 336E-4 

24 

COND = 0.11614 

25 

BETA = 0.01932 

26 

AKN = 18.2 

27 

FRIC = .075 

28 

GO TO 999 

29 

200 CONTINUE 

30 

. . . THE LUBRICANT IS 

31 

VIS1 = 12.76 

32 

VIS2 = 3.2 

33 

A = 10.214494 

34 

B = 3.698398 

35 

RH060 = .9526 

36 

G = 7.092E-4 

37 

COND = 0.15218 

38 

BETA = 0.013168 

39 

AKN = 18.2 

40 

FRIC = .045 

41 

GO TO 999 

42 

300 CONTINUE 

43 

• ■ . THE LUBRICANT IS f 

44 

VIS1 = 25.4 

45 

VI S2 = 4.13 

46 

A = 11.451954 

47 

B = 4.112963 

48 

RH060 = 1.2006 

49 

G = 7 . 47 E-4 

50 

COND = 0.11924 

51 

BETA = 0.016798 

52 

AKN =24.9 


53 


FRIC = .070 

54 


GO TO 999 

55 

400 

CONTINUE 

56 


THE LUBRICANT IS A MIL-L-23699 

57 


VIS1 = 28. 

58 


VIS2 = 5.1 

59 


A = 10.207208 

60 


B = 3.655059 

61 


RH060 = 1.0102 

62 


G = 7.952E-9 

63 


COND = 0.15218 

64 


BETA = 0.016089 

65 


AKN = 18.2 

66 


FRIC = .070 

67 


GO TO 999 

68 C 


LUBRICANT ADDED FOR CRYOGENIC USE 

69 C 


E.S. ARMSTRONG, LEWIS RESEARCH CENTER, 

70 

500 

CONTINUE 

71 C 

• » • 

.THE LUBICANT IS SPECIAL El 

72 


V1S1 = .9370 

73 


VIS2 = .9106 

74 


A = 17 . 2097 

75 


B = 6.9177 

76 


RH060 = 1.569 

77 


G = .002826 

78 


COND = .06921 

79 


BETA = .0180 

80 


AKN = 20. 

81 


FRIC = .070 

82 


GOTO 999 

83 C 


LUBRICANT ADDED - LOW VISCOSITY FUEL 

84 C 


C. M. WOODS , LEWIS RESEARCH CENTER , 

85 

600 CONTINUE 

86 C 


. .THE LUBRICANT IS RP-1 

87 


VIS1 = 1.571 

88 


VIS2 = .7592 

89 


A = 13.98969 

90 


B = 5.26787 

91 


RH060 = 0.8065 

92 


G = 7 . 2E-09 

93 


COND = 0.13806 

94 


BETA = 0.007278 

95 


AKN = 20. 

96 


FRIC = .070 

97 

999 CONTINUE 

98 


CONVERT TO ENGLISH INITS IF REQUIRED 


<♦- 10-86 


10/28/89 



99 

100 

101 

102 

103 

104 


I Ft IMET .EQ. 1 ) GO TO 1000 

RH060 = RH060 * 0.036127 

G = G x 0.5555555 

VIS1 = VIS1 X 1.55E-3 

VIS2 = VIS2 x 1.55E-3 

COND = COND x 0.5777 


105 1000 CONTINUE 

106 RETURN 

107 END 
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Appendix B 
Symbols 


A.B 

constants in Walther's equation 

A, 

generic expression ior /f| and ^2 

A\ 

parameters in the viscoelastic relationship 

a 

semi major length of Hertzian contact 
ellipse, in. 


exponents of the nondimensional traction 
equation 

B, 

lubricant-dependent constant in Barus 
model, °C/Pa 

b 

semi mi nor width of Hertzian contact 
ellipse, in. 

b t ,b. 

coefficients in geometric fit 

bub, 

coefficients in exponential fit 

C 

constants in the Fresco model 

c„, 

general constant of proportionality 

c„ 

constant in the viscoelastic relationship 

C q 

constant in the film thickness relationship 

C,.C\ 

constants of proportionality in traction 
equations 

c,,c 4 

constants in traction equation 

<>■ 

specific heat of fluid, Btu/°F-in. 3 

D„ 

pressure solidification temperature in 
Barus model, °C 

d[)M\ 

coefficients in polynomial fit 

E’ 

effective elastic modulus, lb/in.“ 

G 

fluid coefficient of thermal expansion, l/°< 

h 

minimum film thickness, in. 

k,k; 

conductivity of fluid film, ft-lb/ft F-scc, 
W/m °C, respectively 

»i 

constant in film thickness relationship 

p.p’ 

pressure, lb/in. 2 ; Pa, respectively 

p„ 

maximum Hertzian contact pressure, Ib/in 

AP..\ 

function of P () in viscoelastic relationship 


parameter in the viscoelastic relationship. 
lb/in. 2 

*.v 

curvature sum in .v-direction (rolling 
direction), in. 

T 

fluid temperature, °F 

Tr 

experimental traction force, lb 

Tr 

average experimental traction force per 
unit length, lb/in. 

u 

speed viscosity product 


U ]t U 2 

rolling velocity of disk 1 and disk 2, 
respectively, in. /sec 

w. 

sliding velocity, in. /sec 

u* 

sliding velocity at which ja* occurs, in. /sec 

U s /Us 

nondimensionalized sliding velocity 

V 

roiling velocity (entrainment velocity), in. /sec 

Vo 

constant in the viscoelastic relationship, in. /sec 

Xb.Yr 

shape parameters for the predicted traction 
curve 

X 

coordinate in contact plane in direction of 
rolling, in. 

X * 

coordinate of nondimensional curve for u s /u* 

y 

coordinate in contact plane normal to 


rolling, in. 

Y + 

coordinate of nondimensional curve for m-V* 

z 

coordinate normal to contact plane, in. 

a 

pressure-viscosity coefficient, in. “/lb 

(3 

temperature-viscosity coefficient, 1/ R 

0 ' 

secondary pressure-viscosity coefficient, 
in. 2 /lb 

73 

exponent in the relationship of fi* as a 
function of V 

rjihVi) 

ambient absolute viscosity at inlet tempera- 
ture, lb-sec/in.“; cP, respectively 

V<WllA 

effective absolute viscosity in contact, 
lb-scc/in. 2 ; Pa-sec, respectively 

v\e\ 

ambient absolute viscosity at inlet 
temperature. Pa-sec 

0 

fluid inlet temperature, °C 

X() 

parameter in the viscoelastic relationship 

E 

traction coefficient (tractive force/normal load) 


nondimensionalized traction coefficient term 


maximum traction coefficient 

(91’)* 

maximum nondimensionalized traction 
coefficient term 

V 

ambient kinematic viscosity at inlet 
temperature, eSt 

P*Pf 

fluid density, g/cnv; lb/in. \ respectively 

0 

composite surface roughness, in. 

T.x 

fluid shear stress in rolling direction, lb/in. - 

('/'1I" 2 ) 

nondimensionalized sliding velocity term 

do 1 )* 

nondimensionalized sliding velocity term at 


(*iT 
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